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INTRODUCTION 

The deconfinement phase transition in dense matter, a manifestation of asymptotic 
freedom of quantum chromodynamics (QCD), may lead to stable quark matter cores 
in massive compact stars. The deconfined ("liberated") quarks occupy continuum of 
scattering states, which in the low-temperature and high-density regime fill in the modes 
in a Fermi sphere. The behavior of quark matter in this regime resembles that of 
less exotic low-temperature systems found in condensed matter {e.g. electron gas or 
ultracold atomic vapor) or hadronic physics {e.g. nuclear or neutron matter) [1]. In 
analogy to these system the attractive interaction between quarks, mediated by the 
gluon exchange, which is responsible for the bound state spectrum of QCD, leads to 
quark superconductivity and superfluidity. This phenomenon is commonly refered to 
as color superconductivity [2, 3]. Interestingly, the inhomogeneous superconducting 
phases originally proposed in condensed matter and nucleonic system find their natural 
realization in compact stars when the /3 -equilibrium condition is imposed [4]. 

In the following, I will review and summarize some recent results concerning the as- 
trophysical implications of the hypothesis of deconfinement and color superconductivity 
in quark matter. This contribution is not intended to be a comprehensive review of the 
field and is centered mainly upon my recent research. 



EQUILIBRIUM AND STABILITY 

While we are confident that neutron stars, which are observed as isolated pulsars or 
X-ray sources in binaries, contain nuclear matter in their interiors, the existence of 
hybrid stars with central cores made of quark matter is not obvious. Qualitatively, the 
uncertainty arises because deconfined quark matter has larger number of degrees of 
freedom (color) than the ordinary nuclear matter and, consequently, its formation leads 



to a softening of the equation of states. Such softening may be fatal to the gravitational 
stability of the whole system; in fact, in many models of quark matter equation of state 
no stable configurations are found - the stars collapse into a black hole. A less severe 
uncertainty arises from the unknown physics of deconfinement and the density at which 
it occurs. The central densities of massive neutron stars are at about five times the nuclear 
saturation density or more. The nucleons will hardly remain intact at such high densities 
and the most plausible outcome of compression is a quark-gluon-fluid. 

As an illustration, I start with the study of ref. [5], where the problem of equilibrium 
and stability of non-rotating and rapidly rotating hybrid stars was addressed. This work 
describes the quark matter in the Nambu-Jona-Lasinio (NJL) model, which is a low- 
energy, non-perturbative approximation to QCD, that is anchored in the low-energy phe- 
nomenology of the hadronic spectrum. While dynamical symmetry breaking, by which 
quarks acquire mass, is incorporated in this model, it lacks confinement. A procedure 
which overcomes this difficulty in an ad hoc manner is matching the quark equation of 
state to some low-density equation of state of hadronic matter. This can be done either 
via Gibbs or Maxwell construction; both constructions lead to quantitatively similar re- 
sult, therefore we have adopted the latter. In that case, at the deconfinement phase transi- 
tion there is a jump in the density at constant pressure as illustrated in Fig. 1, left panel. 
The low-density nuclear equations of state are based on the Dirac-Bruckner-Hartree- 
Fock approach [6, 7]. The choice of a particular nuclear equation of state is dictated by 
the need to match it to the quark equation of state. It turns out that it should be rather hard 
in order to fulfill this requirement Alternatively, one may consider the density at which 
the deconfinement takes place as a free parameter by adding an effective bag constant to 
the quark matter equation of state. In these, more phenomenological approaches, softer 
nuclear equations of state can be accommodated [8]. If vector interactions are included 
in the NJL model Lagrangian the stability of the models is improved and softer nuclear 
equations of state can be accommodated as well [9]. 

The spherically symmetric solutions of Einstein's equations for self-gravitating fluids 
are given by the well-known Tolman-Oppenheimer-Volkoff equations [6, 7]. A generic 
feature of these solutions is the existence of a maximum mass for any equation of state; 
as the central density is increased beyond the value corresponding to the maximum mass, 
the stars become unstable towards collapse to a black hole. A criterion for the stability 
of a sequence of configurations is the requirement that the derivative dM/dp c should be 
positive (the mass should be an increasing function of the central density). At the point 
of instability the fundamental (pulsation) modes become unstable. If stability is regained 
at higher central densities, the modes by which the stars become unstable towards the 
eventual collapse belong to higher-order harmonics. 

For configurations constructed from a purely nuclear equation of state the stable 
sequence extends up to a maximum mass of the order 2 M Q (Fig. 1, right panel), 
this large value being a consequence of hardness of the equation of state. The hybrid 
configurations branch off from the nuclear configurations when the central density 
reaches that of the deconfinement phase transition. It is seen that a stable branch of 
hybrid stars emerges in the range of central densities 1.3 < p c < 2.5 x 10 15 g cm~ 3 . 
The masses of the quark cores cover the range < M COK /M Q < 0.75 — 0.88 for central 
densities 1.3 x 10 15 < p c < 2 x 10 15 . Thus, the quark core mass ranges from one third 
to about the half of the total stellar mass. 
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FIGURE 1. Left panel. Number density versus pressure for three models. For the models A and Al 
the nuclear (low density) equation of state is the same; for the models A and B the quark (high density) 
equation of state is the same. Right panel. Dependence of the total stellar mass and the mass of the quark 
core in units of solar mass M e on the central density for non-rotating configurations. The lower set of 
curves represents the masses of the CCS quark cores, the upper set - the total masses of the configurations. 
The maximal masses are marked with boxes. 



Rapidly rotating stars 

Millisecond neutron stars can rotate at frequencies which are close to the limiting 
orbital Keplerian frequency at which mass shedding from the equatorial plane starts. 
The Keplerian frequency sets an upper limit on the rotation frequency, since other (less 
certain) mechanisms, such as secular instabilities, could impose lower limits on the 
rotation frequency. The mass versus central density dependence of compact stars rotating 
at the Keplerian frequency is similar to that for non-rotating stars with the scales for 
mass shifted to larger values [5]. The increase in the maximum mass for stable hybrid 
configurations (in solar mass units) is 2.052 — > 2.462 for the model A, 2.017 — > 2.428 
and 2.4174 for the model Al (there are two maxima) and 1.981 — )■ 2.35 for the model 
B[5]. 

The recent discovery of the double-pulsar system PSR J0737-3039 [10, 11] offers a 
unique opportunity to place further bounds on the gross parameters of compact stars by 
a measurement of the moment of inertia of star A, since the spin frequencies and the 
masses of both pulsars are accurately measured. Timing measurements over a period of 
years could provide information on spin-orbit coupling which could be revealed through 
an extra advancement of the periastron of the orbit above the standard post-Newtonian 
advance or in the precession of the orbital plane about the direction of the total angular 
momentum [12]. The dependence of the moment of inertia, /, of configurations on their 
mass (in the case of rotation at the Keplerian frequency) is shown in Fig. 2. Since 
the moment of inertia is independent of the rotation frequency, we have chosen to 
extract it for configurations rotating at the limiting frequency; note that the masses of 
configurations should be rescaled appropriately, if one is interested in the I(M) function 




FIGURE 2. Left panel. Dependence of Keplerian (maximal) frequency of rotation for hybrid stars on 
their mass (in solar mass units). Right panel. Dependence of moment of inertia of the hybrid configurations 
on the mass of configuration (in solar mass units) rotating at the limiting frequency. The conventions are 
the same as in Fig. 1. 



for slowly rotating configurations. It is seen that while an accurate measurement of the 
moment of inertia of pulsar A can discriminate between the two nuclear equations of 
state, it will not be useful for accessing the properties of hybrid stars since the measured 
masses of pulsars are too low: M/M Q = 1.337 for pulsar A and M/M Q = 1.250 for 
pulsar B; as apparent from Fig. 2 observations of heavier objects are needed to obtain 
useful bounds on the moment of inertia of hybrid configurations. The differences ~ 20% 
in the moments of inertia of purely nuclear and hybrid configurations of the same mass 
(in the case of the models A and Al) are within the accuracy that can be achieved in 
measurements similar to those carried out for PSR J0737-3039. 



Gravitational radiation 

It is interesting that direct observational imprint of the size and mass distribution of 
the quark core will be carried by continuous gravitational waves if the core is crystalline- 
solid [13, 14, 15] or superconducting-solid [16, 17, 18, 19]. Indeed, gravitational waves 
emitted by non-axi symmetric rotating compact stars are expected to be in the band- 
width of current gravitational wave interferometric detectors. Gravity waves arise from 
time-dependent quadrupole deformations of masses. In hybrid stars the axial symme- 
try of the star's core can be broken by the solid deformations in CCS matter. Sev- 
eral authors estimated the strength of gravitational wave emission from crystalline- 
color- superconducting (CCS) matter [17, 18, 19]. The characteristic strain amplitude 
of gravitational waves emitted by a triaxial star rotating about its principal axis is 
ho = \6% Gc el^V r , where v is the star's rotation frequency, r is the distance to 
the observer, e = (I xx — I yy )/I zz is the equatorial ellipticity, is the tensor of moment 
of inertia, G is the gravitational constant, c is the speed of light. The elastic deforma- 



tions are assumed to be small perturbation on the background equilibrium of the star. 
The current upper limit for the Crab pulsar /i < 3.4 x 10~ 25 , which is rotating at the 
frequency v = 29.6 Hz at the distance 2 kpc [20]. The key unknown parameter of the 
theory is the breaking strain of CCS core cr max ~ 10 ^ — 10 ^ and the shear modulus, 
which can be computed approximately and scales as /i s hear ~ A 2 , where A is the pairing 
gap in CCS matter [16]. It follows then that ho A . The strain of gravitational 

wave emission turns out to be close to the upper limit quoted above if A = 50 MeV and 
Omax = 10~ 4 . The Crab pulsar's current limit implies a max A 2 ~ 0.25 MeV 2 , assuming 
maximally strained matter. The evolutionary avenues that may lead to such maximal de- 
formations are not known. Pulsars may simply preserve their initial deformations as they 
cool down and solidify. Additional contribution due to the deformations in the crust of 
the star have not been taken into account in our estimates. However, this contribution is 
likely to relevant to low-mass neutron stars [21]. 

COLOR-MAGNETIC FLUX TUBES 

The information gained from measurements of the gross parameters of compact objects, 
such as the mass or the moment of inertia give us only indirect information on potential 
phases of matter in their interiors. The study of rotational, thermal and magnetic evolu- 
tions of neutron stars provide complementary information, which in many cases is much 
more sensitive to the composition of matter. 

This section discusses the response of dense quark matter to magnetic field following 
Ref. [22]. Our discussion will be focused on the formation and dynamics of flux tubes 
in the two-flavor color superconducting (2SC) phase. The intact gauge symmetry in the 
2SC phase is a linear combination of the ordinary electromagnetic and color symmetries; 
this phenomenon is called "rotated electromagnetism" [23]. The associated gauge field, 
the "<2 photon", is a combination of the original photon and one of the gluons. It is 
massless and propagates freely in 2SC quark matter. The orthogonal combination X is a 
broken gauge generator, and the associated magnetic field has a finite penetration depth. 

Depending on the ratio of the X-flux penetration depth to the coherence length of 
the condensate, the 2SC phase may be type-I or type-II with respect to the X magnetic 
field [24, 25]. A superconductor is of type II if it obeys the condition K = X/t, > 1/ v^2, 
where k is the Ginzburg-Landau (GL) parameter, A is the penetration depth, and ^ is 
the coherence length of the superconductor. For typical values of 2SC superconductor 
K"2SC ~ HA//i g . Therefore, 2SC quark matter will be of type II if the pairing gap is 
sufficiently large, A > H q /16. In quark matter we expect \x q ~ 400 MeV, so this only 
requires the 2SC pairing gap to be greater than about 25 MeV. Similar conclusions were 
reached in refs. [24, 25]. 

Type-II superconductivity leads to formation of flux tubes if the cost of inserting a 
flux tube, which has always a positive self-energy (flux-tension) in a superconductor 
is compensated by the interaction energy of the flux tube with external magnetic field 
intensity. The critical intensity (lower critical field) H c \ for the formation of Abrikosov 
flux tubes containing X-magnetic flux is large, of order 10 17 G. However, we argued 
that, when the quark matter cools into the 2SC phase, the process of domain formation 
and amalgamation is likely to leave some of the X flux trapped in the form of flux tubes. 



The exact configuration and density of such tubes depends on details of the dynamics 
of the phase transition, but the density could be within an order of magnitude of the 
density of conventional flux tubes in proton- superconducting nuclear matter. Note that 
the calculations apply to 2SC quark matter in the temperature range T\$c <T <ti T2SC 
where T^sc is the critical temperature for the formation of the 2SC condensate while 
T\sc is the critical temperature of formation of spin one condensate. 

We shall assume that the flux tubes in the 2SC phase are topologically stable. This is 
not obvious, since they are analogous to "Z-strings" of the standard model [26] which 
have been found to be stable only in a small region of the standard model parameter 
space [27], although the stable region may be enlarged when bound states are taken in 
to account [28]. Because there are also differences between the 2SC phase of QCD and 
the Higgs phase of the standard model, a separate stability calculation will be needed for 
the 2SC case. 

The 2SC phase contains three species of gapless fermions: two quarks ("blue up" and 
"blue down") and the electron. These are expected to dominate its transport properties. 
Strange quarks can be included in the analysis as long as their pairing pattern does 
not break the Q gauge symmetry. Muons may also be present, but, like strange quarks, 
their higher mass gives them a lower Fermi momentum so they make a sub-leading 
contribution to the transport phenomena. 

Light ungapped fermions will be scattered from flux-tubes due to the Aharonov-Bohm 
effect; the cross-section of this process is given by [29] 

do = sin 2 (7r/3) 
d& 27rfcsin 2 (#/2)' 

where /3 = q p /q c , where q p is the charge of the scattering particle. For a flux tube 
that arises as a topological soliton in an Abelian Higgs model, q c is the charge of 
the condensate field whose winding by a phase of 2n characterizes the flux tube; k is 
the momentum in the plane perpendicular to the string, and ■& is the scattering angle. 
Aharonov-Bohm scattering has several important features: (i) The cross-section vanishes 
if /3 is an integer, but is otherwise non-zero, (ii) The cross section is independent of the 
thickness of the flux tube: the scattering is not suppressed in the limit where the symmetry 
breaking energy scale goes to infinity, and the flux tube thickness goes to zero, (iii) The 
cross section diverges both at low energy and for forward scattering. 

The values of j8 for scattering of the fermions that are ungapped in the 2SC phase 
off a flux tube containing magnetic flux is given, to lowest order in the electromagnetic 
coupling constant a, by 

/l a 1 a 1 a 1 a a a\ 

p - dia g^2 + ^'2"2^ , 2"2^ , 2 + 2^'" a? ~ h ~~aJ () 

in the basis defined as i/a = {ru,gd,rd,gu,bu,bd,e~), where "rw" means the red up 
quark, etc. is the electron and a s ~ 1 is the strong coupling constant. We conclude 
that the gapped quarks have /3 close to 0.5, which means that they have near-maximal 
Aharonov-Bohm interactions with an X-flux tube. Among the lighter (and hence more 
phenomenologically relevant) fermions, the Q-neutral bd has zero Aharonov-Bohm 



interaction with the flux tube, while the bu and electron have the same Aharonov- 
Bohm factor sin(^j8 fc ") = sin(7Tj8 e ) m —iz{a/a s ) . The scattering of electrons and blue 
quarks of the flux tubes in the 2SC cores of compact stars contributes to the transport 
coefficients of matter and acts as a dissipative force on the flux tubes. The relaxation 
time for incident light fermion on a flux tube is given by [22] 

T -i = -^sinVA-)- (3) 

PFi 

where n v is the flux density and p Fi is the fermion momentum. It is easy to understand 
the this result. It is of the standard form for classical gases t _1 = cno, where c = 1 is the 
speed of the particles, n = n v is the density of scattering centers, and o °< sin 2 (7r/3)/ pf is 
the cross section for Aharonov-Bohm scattering. One of the blue quarks, the blue down 
quark, has no A-B interaction with the flux tubes (/3 = 0). The other two, the electron and 
blue up quark, have identical A-B factors although their Fermi momenta are different. 

Because the ambient magnetic field in a neutron star is below the lower critical field 
required to force X-flux tubes into 2SC quark matter, the trapped flux tubes will feel a 
boundary force pulling them outwards. The outward force per unit length on the flux tube 
is = r£x I (R 2 — r 2 ) where ex is the energy per length of X flux tubes (flux tension), R 
is the radius of the 2SC core, r is the radial coordinate in cylindrical coordinates. Other 
contributions to the magnetic field from outside of the 2SC core have been neglected. 
Substituting the flux tension we obtain [22] 

fb ^h{w^7 2 )^> K ^ (4) 

where Kx is the Ginzburg-Landau parameter. This force will be balanced by the drag 
force ("mutual friction") on the moving flux tube due to its Aharonov-Bohm interaction 
with the thermal population of gapless quarks and electrons, and also by the Magnus- 
Lorentz force. An estimate of the timescale for the expulsion of X flux tubes from a 2SC 
core on the basis of balance of forces acting on the vortex shows that it is of order 10 10 
years [22]. Therefore, it is safe to assume that the magnetic field will be trapped in the 
2SC core over evolutionary timescales in the absence of other external forces. 

FINAL REMARKS 

There has been substantial progress in understanding the astrophysics of compact ob- 
jects containing quark matter in their interiors over the past years. However much re- 
mains to be done. Further work is needed in developing methods to treat quark matter 
in non-perturbative regime relevant to compact stars. Finding its ground state in the cold 
and dense regime and within realistic models is another long-standing challenge. Astro- 
physics offers several avenues for discerning the state of the matter hidden inside com- 
pact stars. Some of them such as measurements of their gross parameters, gravitational 
radiation and magnetism have been discussed above. Other sensitive tools provided by 
nature are the thermal cooling and anomalies in rotation of compact star. These issues 
are discussed elsewhere [30]. 
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